home *** CD-ROM | disk | FTP | other *** search
/ CD ROM Paradise Collection 4 / CD ROM Paradise Collection 4 1995 Nov.iso / program / nrpas13.zip / RK4.PAS < prev    next >
Pascal/Delphi Source File  |  1991-04-29  |  941b  |  34 lines
  1. PROCEDURE rk4(y,dydx: glnarray; n: integer; x,h: real; VAR yout: glnarray);
  2. (* Programs using routine RK4 must provide a
  3. PROCEDURE derivs(x:real; y:glnarray; VAR dydx:glnarray);
  4. which returns the derivatives dydx at location x, given both x and the
  5. function values y. The calling program must also define the types
  6. TYPE
  7.    glnarray = ARRAY [1..nvar] OF real;
  8. where nvar is the number of variables y. *)
  9. VAR
  10.    i: integer;
  11.    xh,hh,h6: real;
  12.    dym,dyt,yt: glnarray;
  13. BEGIN
  14.    hh := h*0.5;
  15.    h6 := h/6.0;
  16.    xh := x+hh;
  17.    FOR i := 1 TO n DO BEGIN
  18.       yt[i] := y[i]+hh*dydx[i]
  19.    END;
  20.    derivs(xh,yt,dyt);
  21.    FOR i := 1 TO n DO BEGIN
  22.       yt[i] := y[i]+hh*dyt[i]
  23.    END;
  24.    derivs(xh,yt,dym);
  25.    FOR i := 1 TO n DO BEGIN
  26.       yt[i] := y[i]+h*dym[i];
  27.       dym[i] := dyt[i]+dym[i]
  28.    END;
  29.    derivs(x+h,yt,dyt);
  30.    FOR i := 1 TO n DO BEGIN
  31.       yout[i] := y[i]+h6*(dydx[i]+dyt[i]+2.0*dym[i])
  32.    END
  33. END;
  34.